EXISTENCE, UNIQUENESS AND MULTIPLICITY OF 
ROTATING FLUXON WAVES IN ANNULAR JOSEPHSON 

JUNCTIONS 
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Abstract. We prove that the equation modelling an annular Joseph- 
son junction has a rotating fluxon wave solution for all values of the 
parameters. We also obtain results on uniqueness of the rotating fluxon 
wave in some parameter regimes, and on multiplicity of rotating fluxon 
waves in other parameter regimes. 



1. Introduction 
The equation modelling a long Joscphson junction is [7] 

(1) <f>tt + &4>t + sin(0) = <fi xx + f3cj> xxt + 7, 

where a > 0, (3 > 0, 7 > (a, (3 describe dissipative effects, while 7 describes 
the applied bias current). In the case of annular geometry the periodicity 
condition 

(2) cf)(x + L,t) =</)(x,t)+2ir, 

is imposed, where L is the circumference of the junction. 

Definition 1. A rotating fluxon wave is a solution of (p}, © of the form 

(3) </>(x,t) = 9(kx + wt), 
where the function 6(z) satisfies 

(4) 6{z + 2tt) = 6{z) + 2tt Vz e M. 
Clearly the assumption that satisfies © implies that 

(5) k = j-, 

while the frequency u> is an unknown to be found as part of the solution (in 
Proposition [5] we show that necessarily u> > 0). We note that the time-period 
T of rotation of the wave around the annulus is given by T = —. 
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Substituting © into (JTJ) we obtain 

(6) I3uk 2 9>"(z) + (k 2 - lu 2 )9"(z) - au9'{z) - sin(0(z)) +7 = 0, 

so that rotating fluxon waves are given by solutions (lo, 9{z)) of ([6]), with 6{z) 
satisfying (g]). 

Our aim in this paper is to study the questions of existence and unique- 
ness/multiplicity of rotating fluxon waves. We shall obtain existence of a ro- 
tating fluxon wave for all values of the parameters: 

Theorem 2. For any L > 0, a > 0, (3 > 0, 7 > 0, 0), (0) has a rotating 
fluxon wave solution (0), with lu > 0. 

To prove Theorem [5] and our other results, we use methods of nonlinear 
functional analysis, reducing the problem to a fixed-point equation depending 
on a parameter (lu), to be solved by means of Leray-Schauder theory coupled 
with an auxilliary scalar-valued equation, which we analyze. We have used a 
similar technique to prove existence of travelling waves of the discrete damped 
and dc-driven sine-Gordon equation [5]. The details of the analysis differ, and 
in the discrete sine-Gordon case we could not prove existence for all parameter 
values (which reflects the phenomenon of 'pinning' in the discrete case) - while 
here we are able to do so. 

Several researchers have studied long Josephson junctions, and annular 
ones in particular, experimentally, [21 E E] , numerically [IJ [6] and analytically 

[aimiiiz]. 

Note that in the case (3 — (known as the damped dc-driven sine- 
Gordon equation), © reduces to the damped dc-driven pendulum equation 
which is amenable to phase-plane analysis and is thus well understood, and 
the existence of rotating fluxon waves for all a > 0,7 > is known (see [TO] , 
Section 2.2). The situation for (3 > is different. Most existing analytical 
studies assume that either a, (3 and 7 are small, so that the equation is 
treated as a (singular) perturbation of the sine-Gordon equation (the case 
a = 0, (3 = 0, 7 = of ((T|)), or that (3 is small so that the equation is a 
perturbation of the damped dc-driven sine-Gordon equation. An exception is 
[6j, which studies travelling fluxon waves of kink type on an infinite line, that 
is, solutions of (p} of the form ([3]) with ((4]) replaced by 9(+oo) — 9(—oo) = 2ir. 
Existence results for such waves are proved by a geometric analysis of the 
corresponding three-dimensional phase-space, and a very intricate structure 
of solutions is revealed. Here, as we noted, our approach is functional-analytic 
rather than geometric, leading to the general existence result of Theorem [2) 

Our method also allows us to derive some further information about the 
set of rotating fluxon waves. In the physical literature, a common and useful 
way to describe the behavior of a system is the TV characteristic' - in this 
case the 7 — lu (bias current-frequency) characteristic: 
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Definition 3. The bias-frequency characterstic is the set x C (0, oo) x (0, oo) 
of all pairs (a;, 7) for which ([I]),© has a rotating fluxon wave ([3|), with fre- 
quency LO. 

The set x depends, of course, on the parameters L,a,f3. The following 
theorem shows that the set x is 'large'. 

Theorem 4. For any L > 0, a > 0, /3 > 0, the bias-frequency characteristic 
X has a subset xo C X which has the following properties: 

(i) Xo is connected. 

(ii) For any w > there exists 7 > with (w, 7) S xo- 
(Hi) For any 7 > there exists u> > with (w,7) S Xo- 

H (0,0) exs. 

(W X \ Xo * s a bounded set. 

Note that Theorem [5] follows at once from part (hi) of Theorem [U 
We shall also obtain some bounds on the frequency u of rotating fluxon 
waves (see section [5]): 

Proposition 5. For any rotating fluxon wave of we have w > and 

(7) < w < — . 

We note that since we are assured that w > 0, the lower bound in ([7]) is 
nontrivial only when 7 > 1. 

The above results are completely general in that they hold for all values 
of the parameters. We also obtain some further results under some restrictions 
on the parameters. While we do not expect our conditions on the parameters, 
which are imposed by our methods of proof, to be sharp, these more restricted 
results expose some of the phenomena which can occur and point the way to 
further investigation (analytical or numerical) by raising various questions 
which remain unresolved. 

We do not know whether in general we can take Xo — X m Theorem 
IU or in other words whether the characteristic x is connected. In Section [3] 
we shall provide a positive answer to this question, and also obtain further 
information about the bias-frequency characteristic x, for a certain subset of 
the space of parameters L,a,f3: 

Theorem 6. Assume L > 0, a > 0, (3 > 0. If 

(8) k = T >1 

and 

(9) a + flk 2 > V2k 

then the bias-frequency characteristic x can be represented in the form 

(10) x = {(^,7(w))|w>0} 
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where 7 : [0, 00) — > [0, 00) is continuous and real-analytic in (0,oo), and 
satisfies 

(11) 7(0) = 0, 

(12) lim 7(0;) = +00. 

UJ — > + oc 

In particular, \ is connected, and for a given bias 7 > the frequencies of 
rotating fiuxon waves are the solutions lu of the equation 

(13) 7(0; ) = 7. 

Let us note that in the case (3 = 0, part of the above theorem is not 
true: while \ is connected and can be parameterized as in (|10[) (for any k > 0, 
a > 0; see [10], sec. 2.2), the function 7 is not analytic in the case (3 = Q, in 
fact it has a cusp at ui = k (see [10], fig. 2), and the value lu = k separates 
the 'low-voltage' and 'high-voltage' regions. Theorem [6] shows that the term 
4>xxt in ([1]) has a smoothing effect on the curve x ( a t least in the case that 
©,([1]) hold). When we prove Theorem[S]we will point out why the assumption 
(3 > is crucial, which leads to this difference between the damped dc-forced 
sine- Gordon equation and (fTj). 

A natural question is whether the rotating fiuxon wave whose existence 
is guaranteed by Theorem [21 is unique (for given L, a, (3,"f). In Section [4] we 
obtain a uniqueness result in the case that 7 is sufficiently large. 

Theorem 7. Fix any L > 0, a > 0, (3 > 0. Then for sufficiently large 7 > 0, 
there is a unique rotating fiuxon wave. 

On the other hand, in Section [5] we show that uniqueness does not hold 
in general, by proving that for certain parameter ranges there exist at least 
three rotating fiuxon waves with different frequencies. 

Theorem 8. Assume and replace (G|) by the stronger assumption 

(14) (3k > V2. 

Then there exists a value 70 > such that given any < e < ^7o 7 there exists 
«o > so that for any a € (0, ao) and any 7 € (e,7o — e), {!]), 0) has at 
least three rotating fiuxon solutions with distinct frequencies. 

We can expect this multiplicity of rotating fiuxon waves to lead to bista- 
bility, jump and hysteresis phenomena as the bias 7 is varied. We do not 
know whether the result of Theorem [8] is valid for all values of the parameters 
L, a, (3, or whether, on the contrary, there are some values of these parameters 
for which there is a unique rotating fiuxon wave for all 7 > 0. 

We mention here that all our results and proofs in this paper remain 
valid, with minor changes, if the nonlinearity sin(0) in ([T]) is replaced by p(4>), 
where p is an arbitrary C 1 27r-periodic function, satisfying p{9)d9 = 0. 
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An important issue which we do not address here is the question of 
stability. For example: is it true that for any value of a > 0, (3 > 0, 7 > 
there exists at least one rotating fluxon wave of (JTJ) which is at least locally 
asymptotically stable? 



2. Existence of rotating fluxon waves 

In this section we prove Theorem[4j which in particular implies Theorem 
[2] To do so we first recast our problem into a functional-analytic form, so that 
we can apply Leray-Schauder theory. 

Multiplying both sides of ([6]) by 0'(z) and integrating over [0, 2tt], taking 
(HJ) into account, we obtain 

(3k 2 [ {8"{z)) 2 dz + a ( (9'{z)) 2 dz 



27T7- 

By our assumption that a > 0, (3 > 0,7 > this implies that ui > 0, as stated 
in Proposition [5l and thus we may set 

(15) X=- 

U) 

We also set 

(16) 6{z) = z + u(z), 
with u(z) satisfying 

(17) u{z + 2tt) = u(z) Mz e M, 
so that we can rewrite © as 

(18) Pk 2 u"'(z) + (Afc 2 - j)u"(z) - au'(z) - Asin(z + u{z)) + A7 - a = 0. 

Rotating fluxon solutions of fl}, ^ thus correspond to pairs (X,u(z)) with 
A > and u(z) satisfying (fT7|) . (fT5|) . 

By integrating (fTS]) over [0, 2ir], taking (| 1 T[) into account, we obtain 

1 / ot 

(19) — J sm(s + u(s))ds = 7- -, 

so that we can rewrite (TT51) as 



(3k 2 u"'{z) + (Afc 2 - j)u"{z) - au'(z) 



1 r 2n 

(20) - Asin(z + u(z)) + A — / sin(s + u(s))ds = 0. 

27r Jo 

We note that if (A, u(z)) is a solution to our problem then so is (A, u(z + c) +c) 
for any constant c (this of course corresponds to the time-invariance of CU), 
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hence by adjusting c we may assume that 

r 2TT 



(21) 



u(s)ds = 0. 



We thus seek pairs (X,u(z)), satisfying (fTf|). (fT9| . ([20]) and ([21]) . 

Multiplying (f20|) by 1 + u'(z) and integrating over [0, 27r] we obtain 

Lemma 9. Any solution (A, u) of {17jj , ifgflp and satisfies 

(22) A / sin(s + u(s))ds = (3k 2 (u"(s)) 2 ds + a (u'(s)) 2 ds. 



In particular, together with (|19]) . ([22)1 implies that 

A7 > a, 

which gives the upper bound of (0. Also, noting that the value of the integral 
on the left-hand side of (fH?)) is at most 1, we have 



(23) 



7-l<p 



which implies the lower bound in ([7]). We have thus proven Proposition [5j 
We denote by X, Y, Z the Banach spaces of real- valued functions 
X = ju G i? 3 [0,27r] m w (0) = u (2) (2tt), < i < 2, ^ u(s)ds = o| , 

y = lu G ^[0,2^] ^ u(s)ds = q|, 



Z= «€ L 2 [0,2tt] 



2tt 



u(s)ds = > , 



with the norms 



/ 1 



2;r 



NU = I (u'"(s)fds 



2 IT 



H|z 



1 

27 



271 



(u(s)) 2 ds 



Assuming v G i? 1 [0, 27r] with Jg 2 ^ v(s)ds = 0, developing u in a Fourier series 
and using the Parseval identity one obtains 



(24) 



2tt 



2tt 



(v(s)Yds< / (v'(s)Yds, 



known as Wirtinger's inequality, implying the following inequalities which we 
will often use: 



Idly < llullx € X, \\u\\ z < \\u\\ Y Vm G y. 
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We define a family of linear mappings L\ : X — > Z (A ^ 0) by 
(25) L A (u) = (3k 2 u'"{z) + (Afc 2 - - cra'(«) 

and a nonlinear mapping N : Z — » Z by 



(26) 



1 



i\T(u) = sin(z + u(z)) — — / sin(s + u(s))ds 



271 



Solving ([IT]). |20|) , (|2T|) is equivalent to solving 

(27) L A (u) = AiV(tt), ueX, 

so that our problem reduces to: find (A, it) € (0, oo) x X such that (|27ll and 
(I19[) hold (we note that, by a simple bootstrap argument, any solution tiel 
of (f2"T)) is in fact C°°). In other words, defining the nonlinear functional $ : 

(0, oo) x Z -> M by 

a 1 /" 27r 

(28) $(A,n) = ~ + — / sin(s + u(s))ds, 

A 27T J 

we need to solve the pair of equations l|27j) and 

(29) *(A,«)=7, 

for the unknowns (A, u) € (0, oo) x X. 

We derive some estimates on L\ which will be useful. Let u € X and 
write it as a Fourier series 

u(z) = V" a/e li2 \ 



Noting that 



L x {e ilz ) =- (la + (3k 2 l 3 )i + l 2 (Xk 2 - A" 1 ) 



ilz 



so that 



L^(e ilz ) = - (la + (3k z l :i )i + l 2 (Xk 2 - A" 1 ) 



s 
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we obtain 

\\LxHu)\W = 



^2 a l [i la + Pk 2 l 3 )i + l 2 (\k 2 - A _1 )J \ llz \\ 



Z#0 



= (]>>| 2 [(a + /3fc 2 ^ 2 ) 2 + ; 2 (Afc 2 - A" 1 )' 



In i 



Z#0 



< [(a + /3fc 2 ) 2 + (Afc 2 -A- 1 ) 2 ] _3 [^ 



so that 
(30) 



(a + /3fc 2 ) 2 + (Afc 2 -A- 1 ) 2 || tt || 



H^lky < (a + /3fc 2 ) 2 + (Afc 2 -A- 1 ) 2 



In particular (dropping the first term on the right-hand side of (|30|) ) 



(31) 



< 



\X 2 k 2 - If 



We can also obtain 



||iA X («)IU = (E^|a/| 2 [(^ + /3fc 2 Z 3 ) 2 + Z 4 (Afc 2 



< 



|A 2 fc 2 - 1| 



A 



z#o 



|A 2 fc 2 - 1| 



-1\2 



so that 
(32) 



A 



|A 2 fc 2 - If 
Defining F : (0, oo) x Z -> Z by 

(33) F(X,u) = XL' 1 o N(u), 
we may rewrite (|27p as 

(34) u = F(X,u), ueZ. 

We note that since L^ 1 maps Z onto X, any solution ugZof (f3"4")l is in fact 
in X. 

Since || sin(z + u(z))\\ L 2 < 1 and N(u) is the orthogonal projection of 
sin(z + u(z)) £ L 2 [0,2tt] into Z, we have 



(35) 



\\N(u)\\ z <l VueZ. 
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Using ((35]) and ([30]) . we have the bound: 
(36) 

\\F(X,u)\\z < \\F(X,u)\\ Y < A[(a + /3fc 2 ) 2 + (Afc 2 -A- 1 ) 2 



Vm € Z, A > 0. 



We extend the definition of F(X,u) by setting F(0,u) = for all u 6 Z 
(this is necessary because La is undefined for A = 0). The bound ([56]) shows 
that the extension is continuous, and moreover that the mapping F takes 
bounded sets in [0, oo) x Z into bounded sets in Y , hence (since the embedding 
of Y into Z is compact) into compact sets in Z, so that F : [0, oo) x Z — > Z 
is a compact mapping. 

We define 

£ = {(A, it) G (0,oo) x Z | u = F(A,u)}. 
The solutions of (f2T)) . (f!?9")) , are the solutions of the equation 

(37) *(A,«)=7, (A,u)GS. 

In other terms, recalling (|15|) . the bias- frequency characteristic x defined in 
the introduction is given by 

(38) x={(",H"-\u)) | (u>-\u)€Z}. 

We will now apply the following version of the Leray-Schauder continu- 
ation principle (see, e.g., [TT], Theorem 14. D): 

Theorem 10. Let Z be a Banach space and let F : [0, oo) x Z — > Z be a 

continuous and compact mapping, with 

F(0,u) = VueZ. 

Define Sc(0,oo)xZ by 

£ = {(A, it) G (0, oo) x Z | u = F(X, u)}. 

Then there exists an unbounded connected set Eq C E with (0,0) £ Eo- 

Since the mapping F that we have defined satisfies the hypotheses of 
Theorem [TUl we have a set Eo as in the Theorem. 
We define xo C x by 

(39) xo - {(«, Hu>-\ v))\u> 0, {lo-\u) e E }, 

and we will show that xo has all the properties claimed in Theorem |U 

Since $ is continuous and Eo is connected, xo is connected, so we have 

(i) of Theorem gj 

Part (ii) of Theorem 2] follows from the following 

Lemma 11. Let Eo be as in Theorem ] 1 (A For any A > there exists u G Z 
such that (A, u) G Eo. 
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PROOF: Since from (l34j) and ([36]) we have 

(40) (A , u)eE ^ N | z < A _^_. 

we conclude that the unboundedness of E is in the A-direction, that is, for 
any M > there exists (A, it) € So with A > M. In other words, denoting by 
A c (0, oo) the projection of E to the A-axis, 

A = {A > | 3u e Z such that (A, u) G E }, 

we have that A contains arbitrarily large values. Since (0,0) G Eo, A also 
contains arbitrary small positive values. Since So, hence A, is connected, we 
conclude that A = (0,oo). □ 

To prove that xo satisfies property (iii) of Theorem |4] (and hence also 
Theorem [2]) , we need to show that 

(41) (0,oo)c*(E ). 

Since So is connected and $ is continuous, ^(So) is also connected, so that 
to prove (|4"Tj) it suffices to show that $ takes arbitrarily large and arbitrarily 
small positive values on So. Since the integral in the definition (|28|) of $ is 
bounded by 1, it is immediate that 

(42) lim $(A, u) = +oo, 

A^0+, (A,u)6£ 

so that $ takes arbitrarily large values on So. We shall prove 
Lemma 12. 

(43) lim $(A,u) = 0. 

A^+oc, (A,ti)eS 

This implies that $ takes arbitrarily small positive values on So, com- 
pleting the proof of (|4"Tj) , and thus of part (iii) of Theorem 0] holds. 

PROOF OF Lemma [TJJ Applying Wirtinger's inequality (|2"4"|) with v = u', 
Lemma |2] implies 
(44) 



f 



2tt 



f 



(A,u)eS A— / sm(s + u{s))ds < {(3k 2 + a)— / {u'{s)) 2 ds 



2n 



2 TV 



From (|34l) and (|31|) we obtain 

A 2 



(45) (A,u)eE =► \\u\\ Y < , yk , 

(|44| together with (|45| implies 



f 



(A, u)sS =^> — / sin(s + u(s))ds 



2 71 



^ (/3fc 2 + a)A' i 
" (A 2 fc 2 -1) 2 
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which implies that 

(46) lim — / sin(s + u(s))ds = 0, 

hence we have (|4U)). □ 

To prove claim (iv) of Theorem IU let u n = — (n > 1), and note that by 
part (ii) of the Theorem, there exists, for each n, 7„ > such that (u> n ,j n ) G 
Xo- By the definition of xo there exists u n € So such that j n — <&{uj~ 1 ,u n ). 
But by Lcmma [T2l we have 7„ — > as n — > oo. Hence we have (w„, j n ) — > (0,0), 
so we have (0, 0) € Xo- 

To prove part (v) of Theorem [U we will use the following two lemmas, 
which will also be used in the following sections. 



Lemma 13. If 



(47) A (a + /3k 2 ) 2 + (Afc 2 - A" 1 ) 2 



< 1 



then F(X, .) : Z —>Y (hence also F(X, .) : Z — > Z) is a contraction. 
PROOF: Computing the Frechet derivative of TV : Z — > Z we have 



N'(u)(v) = cos(z + u(z))v(z) — — / cos(s + M(s))w(s)ds, 

2tt Jo 

so that AT' (it) is the mapping from Z to L 2 [0, 27r] given by w(z) — > cos(z + 
m(z))w(z) followed by an orthogonal projection into the subspace and this 
implies 

(48) ||W / (u)|| ZiZ < 1 VueZ. 



Vu e z. 



From (|30]l. (|33| and gS]) we conclude that 
(49) 

\\D u F{\,u)\\z,y - A||LJ- 1 JV / (u)|| z ,y < A[(a+/3fc 2 ) 2 +(Afc 2 -A- 1 ) 2 
which, by (|47| . completes the proof. □ 

Lemma 14. Fixing a > 0, /3 > 0, i/iere exists Ao > sitc/i i/tai /or any 
A G [0,Ao) i/te equation {3J$ has a unique solution, which we denote by u\. 
The mapping A — > u\ from [0, Ao) into Z is continuous and is real- analytic 
for A G (0, A ). 



PROOF: It is easy to see that we can find Ao > such that (|47|) holds for 
A G [0, Ao), and thus by Lemma [T51 for any A G [0, Ao) the mapping F(X, .) : 
Z — ► Z is a contraction, so that, by Banach's fixed-point Theorem, equation 
(f54"| has a unique solution, which we denote by u\. Real-analyticity of the 
mapping A — > it a follows from the real-analytic implicit function Theorem. □ 
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Returning to the proof of part (v) of Theorem [4] we note that since, 
by Lemma 111) for each A there exists u such that (A, u) G So, and since the 
uniqueness in Lemma [TJ] implies that when A < Ao such a u is unique in E, 
we must have 

(50) < A < Ao, (A, u)eE => (A, u) G E . 
Thus if we set w = ^ we have by ([50]), ([38]), ([39]) 

(51) u > ujq, (w,7)Gx => (w,7)exo- 
Setting 70 = aujQ + 1 and using Proposition [5] and (|5ip we obtain 

(52) 7 > 70, (w,7)£x w > ^o, (w,7) £X =>■ (^,7) G Xo, 
and I[51]) .([g2 ]l mean that 

7) € x \ xo w < cj , 7 < 7o, 
so that x \ Xo is bounded. 

We note for later use that ([22]) implies that 
Lemma 15. 

1 /- 2t 

(A,m)gE, A>0^— / sin(s + u(s))ds > 0. 
27r Jo 

The strictness in the above inequality follows from the fact that, by 
(f2"2"]l . if we had equality this would imply u = 0, but it is easy to check that 
(A, 0) g E for A > 0. 

3. More on the bias-frequency characteristic 

In this section we prove Theorem [5] which shows that for a certain open 
subset of the space of parameters L, a, [3, the bias-frequency characteristic x 
is a connected, real analytic curve, and can be represented as the graph of a 
function: 7 = 7~(w). 

Lemma 16. Assume that a > 0, f3 > and i/ie inequalities and (0) /10/d. 
T/ien for all A > 0, the equation Jff^| ) /ias a unique solution, which we denote 
by u\, and therefore 

(53) E = E = {(A, MA ) I A>0}. 

The mapping A — > ma /rom [0, 00) into Y is continuous, and real-analytic on 
(0,oo). 

PROOF: We use Banach's fixed point Theorem. By Lemma [TBI to show that 
F(X, .) is a contraction from Z to itself for all A > it suffices to show that 
(|47| holds for all A > 0. After some manipulations, and setting fi = A 2 , we 
see that this inequality is equivalent to 

(54) p(n) = (k 4 - l)n 2 + [(a + /3k 2 ) 2 - 2k 2 ]n + 1 > > 0. 
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The condition k > 1 ensures that the quadratic function p(p) has a minimum, 
and we have p(0) = 1 > 0. Moreover the condition {9j ensures that p'(0) > 0, 
which implies p is increasing on [0, oo), so that (154"]) holds. We thus obtain 
a unique solution u\ to (|34[) . The fact that it is real-analytic in A G (0, oo) 
follows from the analytic implicit function Theorem and the fact that the 
function A — > L^ 1 from (0, oo) to the space B(Z) of bounded linear operators 
from Z to itself is real-analytic. □ 

Let us note here that it is precisely the last statement in the above proof 
- the fact that A — > L^ 1 is real-analytic, that breaks down when (3 = 0, and 
is responsible for the fact that in that case (i.e. the case of the damped, de- 
forced sine-Gordon equation) the curve x is not smooth (see [10J . fig. 2). The 
reason that, when (3 — 0, A — > L^ 1 is not real- analytic can be gleaned by 
recalling the definition 1(25]): if j3 = then at A = r the coefficient of the 
highest derivative in the operator L\ vanishes. 

Thus, under the assumptions of Lemma [TBI we have 

X = Xo = fa' 1 )) |^>0}, 

where \1/ : (0, oo) — > K is defined by 

(55) ¥(A) = #(A,u A ), 

so defining 7(0;) = ^(w -1 ) we have ((TP]) . The fact that 7(0;) > for w > 
follows from Lemma [131 (JTTJ) follows from Lemma [121 and (JT2J) follows from 
(]4"2)) . This concludes the proof of Theorem (HI 

4. UNIQUENESS OF ROTATING FLUXON WAVES FOR LARGE BIAS CURRENT 

We now prove Theorem [7] By Lemma 1141 defining Ao as in Lemma 
dH for any A € [0, Ao) the equation ((34]) has a unique solution, which we 
denote by U\. We define the function ^> : (0, Ao) — > M by ((55]) . where $ is the 
functional defined by ((25]) . so that solutions of ((37]) with A € (0, Ao) correspond 
to solutions of 

(56) *(A)=7. 

We show below that there exists some Ai € (0, Ao) such that 

(57) Ae(0,Ai) =*> *'(A)<0. 

We note that if 7 is sufficiently large then, by (O, we have A S (0, Ai) for any 
solution of |56]) . Therefore ((57]) implies uniqueness of solutions of ((56]) , hence 
uniqueness of the rotating fluxon wave, for sufficiently large 7, which is the 
content of Theorem 
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To prove claim ([57]) we note that, since 

/■2jt 



a 1 r 

*(A) = - + — J sin(s + u x (s))ds, 



*'( A ) = ~T2 + / C0S ( S + MA(s))-DAUA(s)ds, 



we have 

a 1 

and using Cauchy's inequality 

(58) *'(\)<-^ + \\D x u x \\ z , 

so to prove (|57|) it suffices to show that ||-Dau a ||z = °( j?) as A — > 0+. In fact 
we shall prove a much stronger statement: 

Lemma 17. We have 

(59) PauaIIz = O(A) asA^0 + . 
PROOF: We differentiate the identity 

L x {u x ) = \N(u\) 

with respect to A, obtaining 

(60) (D\L\)(u\) + L x {D x u x ) = AT(u A ) + AJV'(u A )(£> A u A ). 

We note that, because we assume A <E (0, Ao), where Ao is defined as in 
Lemma [Til AL^ 1 N' (u x ) is a contraction from Z to itself so the inverse of 
I - XL^N'iux) is well defined as a mapping from Z to itself. Moreover we 
have, using (|49j) . 
(61) 

|| [7 - XLfN'(u x )]- l \\ z , z < ... * , Ml = aS A ° + ' 

1 - A||i A N'(u x )\\z,z 

We can rewrite (pD]l in the form 

J X N \".\>! 



D x u x = [I- XL^N'iu^L-^N^) - (D x L x )(u x )], 



which implies 
(62) \\D x u x \\ z 

< ||[/-AL A 1 7V'( UA )]- 1 || z , z [||i A 1 7V( UA )|| z + ||L- 1 (/J A i A )( UA )|| 

Using P5|) we obtain 

H^MHy = / c os 2 (s + »a(*))(1 + <{s)) 2 d S ) 5 



- (^/* r ( 1 + «A(*)) 2 d») i = [l + IMI^]'< [l+( 



A 2 



A 2 fc 2 
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so that 

(63) \\N(u x )\\y = 0(1) as A ->() + . 

From ((32)) and (|63|) we have 
(64) 

||£^JV(u A )|U< WL^NiuxjWx < \\Li 1 \\ Yt x\\N(ux)\\ Y = 0{\) as\^0 + . 
From the definition (|25|) of L\ we have, for all u E X 

(D x L x )(u)=(k 2 + ±)u"(z), 

hence 

(65) || (Ma) (u) || y - (k 2 + ^)\\u\\x V«el. 
From ([64]) we have 

(66) ||«a||jc =M\Lx 1 N(u x )\\x = 0(A 2 ) asA^0+, 
so using (|65|) . (|66|) 

(67) ||(D A L A )(u A )||y=0(l) as A -> 0+, 
and together with (f3"2")) 

(68) ||£^(D A L A )(u A )|| z < \\Lf(D x L x )(u x )\\ x = 0(A) as A -> + . 
Using ([6T ]> . ([62 ]) . (|64 |) and (p] l we obtain (59]). □ 

5. MULTIPLICITY OF ROTATING FLUXON WAVES 

In this section we prove Theorem [8] which shows that for some values of 
the parameters L,f3, if a > is sufficiently small then there exists an interval 
of 7's for which there exist at least three rotating fluxon waves. 

Lemma 18. Assume that (3 > and |J^[ ) hold. Then for any a > 0, 
A > 0, the equation \3J$ has a unique solution, which we denote by u atX . 
The mapping (a, A) — > u a X from [0, oo) x [0, oo) into Z is continuous. 

As in the proof of Lemma HH this follows from the Banach fixed point 
Theorem applied to Note that (HI implies that © holds for all a > 0. 

Since in the following argument the dependence on a is important, we 
display it in the notation (of course u a ^ x also depends on the parameters k and 
(3, but since we shall regard these as fixed we can suppress this dependence). 

PROOF OF Theorem [8j Under the conditions of Lemma [18] we can define 
* a (A) as in (55]) . now defined for all A <E [0, oo), so that rotating fluxon waves 
correspond to solutions of 



(69) 



* tt (A)=7, A>0. 
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We write * Q (A) in the form (see [[28])) 

(70) * a (\) = j + r a (\), 
where 

1 f 27r 

= / sin(s + u a \(s))ds. 
27T Jo 

We note some properties of the function T Q (A) (valid for any a > 0). 
Since it Qj o = we have 

(71) T Q (0) = 0. 
From (|4"6")l we obtain: 

(72) lim T Q (A) = 0. 
Lemma [T5] implies 

(73) T a (A) > VA > 0. 

(|7T|) . ([72^1 . (173^) , imply that T Q (A) attains a positive global maximum at some 
A > 0. In particular this is true for a = 0, and we denote the global maximum 
of To(A) by 70, and assume that it is attained at Ao > 0: 

(74) T (Ao)=7o = niaxT (A). 

We now fix an arbitrary < e < ^70. By (f7"Tj) . ([72|) we can choose < Ai < A 
and A2 > Ao so that 

(75) T (A 1 )<|, T (A 2 )<^. 

Since the dependence of T Q (A) on a is continuous, we can choose ao > 
sufficiently small so that 

(76) |T„(A) - T (A)| < I Vae [0,a o ],Ae [0,A 2 ]. 
From dT3I) ,dZnD and J7HD we obtain, for all a £ [0, a ], 

(77) T Q (A )>7o-f, 

(78) T Q (A 1 )<^e, T Q (A 2 )<Je. 
Moreover we can choose «o > sufficiently small so that 

From ([70]) and (|79|) we then obtain 

|# a (A) - T Q (A)| < I Vae [0,a o ],A> Ax. 
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Together with f77| l .([78 |l this gives, for all < a < a , 
(80) * Q (A ) > 7o -e, 



(81) 



* a (Ai) < e, 




lim \& a (A) = +oo. 



*a(A 2 ) < e. 



Thus by continuity of \& Q on (0, oo), if we take 7 <E (e, 70 — e) then by and 
([81]) the equation (|69j) has a solution AG (0, Ai), by (j8l"j) and ([80]) the equation 
([69]) has a solution A £ (Ai, Ao), and by (j80|) and ([82]) the equation (|69|) has a 
solution A G (Ai, A 2 ). We have thus found, under the stated assumptions, three 
solutions of (f69|) . corresponding to three rotating fluxon waves with different 
frequencies. □ 
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